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ABSTRACT 

An  experimental  test  program  was  undertaken  to  characterise  the  inelastic  response  of 
D6ac  steel  under  differing  cyclic  loading.  A  unified  constitutive  model  was  used  to 
describe  material  behaviour.  Constants  for  the  model  were  extracted  from  the  test 
results  and  numerical  model  predictions  were  compared  with  the  experimental  results. 
The  tests  revealed  that  D6ac  steel  saturates  after  one  cyclic  loop  and  a  change  in 
modulus  was  observed  after  rmloading  from  a  tensUe  load.  Also  described  are 
improvements  in  the  testing  of  uniaxial  specimens,  focusing  on  techniques  for 
specimen  alignment  to  minimise  out  of  plane  bending  effects. 
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Characterisation  of  D6AC  Steel  Using  A 
Unified  Constitutive  Model 


Executive  Summary 


Every  7  years  the  F-111  Aircraft  undergoes  a  Cold  Proof  Load  Test  (CPLT)  as  part  of 
the  Structural  Integrity  Program.  During  this  loading,  local  regions  of  the  Wmg  Pivot 
Fittinff  experience  a  high  level  of  plastic  deformation.  A  consequence  of  particular 
concern  is  that  the  Fuel  Flow  Vent  Hole  Number  13  (FFVH#13)  experiences  faPgue 
cracking  due  to  the  residual  stresses  which  are  left  m  the  structure  after  the  CPIT 
cycle.  The  Royal  AustraHan  Air  Force  (RAAF)  manages  the  structural  integrity 
111  fleet  by  safety  by  inspection,  with  inspection  intervals  derived  from  a  durability 
and  damage  tolerance  analyses  (DADTA).  As  part  of  the  calculation  of  the  inspectaon 
interval  for  FFVH#13,  a  detailed  knowledge  of  the  residual  stress  field  imposed  by  a 
single  CPLT  cycle  is  required.  AMRUs  role  is  to  accurately  determine  the  elastic  and 
cycHc  residual  stresses  as  input  to  this  DADTA. 


Recent  full  scale  wing  tests  conducted  at  AMRL  indicate  that  at  peak  load  the  lower 
inboard  comer  of  FFVH#13,  where  instances  of  cracking  have  occurred,  experiences 
strain  levels  three  times  that  of  yield.  To  accurately  analyse  the  effects  of  such  high 
strain  levels,  a  Finite  Element  analysis  of  the  critical  region  was  required.  This  analysis 
must  correctly  represent  the  material  behaviour  under  the  non-symmetric  cyclic  loads. 
Classical  plasticity  models  have  inherent  difficulties  when  dealing  with  cyclic  loading 
and  therefore  alternative  methods  have  been  employed.  One  such  method  used  at 
AMRL  is  a  state  variable  constitutive  material  approach,  originahy  developed  by 
Ramaswamy,  Stouffer  and  Bodner,  and  subsequently  jointly  enhanced  by  AMRL,  the 
University  of  Melbourne  and  Monash  University  over  recent  years. 


This  report  describes  an  experimental  test  program  which  characterises  the  inelastic 
response  of  D6ac  steel  under  differing  cyclic  loading.  Unified  constitutive  material 
model  constants  were  extracted  from  the  test  results  and  refined  using  a  FORTRAN 
program  which  represents  the  uniaxial  test  condition.  These  refined  constants  when 
apphed  to  the  constitutive  model  successfully  predicted  the  behaviour  of  D6ac  steel 
under  different  cyclic  loading  conditions  including  the  CPLT  cycle. 


This  constitutive  model  forms  part  of  the  input  to  the  Finite  Element  solutions  used  in 
the  detailed  plasticity  analyses  of  the  F-111  FFVH#13.  The  elastic  and  residual  stress 
fields  foimd  in  this  analyses  will  then  be  used  in  tihe  DADTA  as  an  important 
ingredient  in  deterrnining  the  inspection  interval  for  this  critical  location. 
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1.  Introduction 


The  Aeronautical  and  Maritime  Research  Laboratory  (AMRL)  has  the  key  role  of 
supporting  the  Royal  Australian  Air  Force  (RAAF)  in  the  structural  integrity 
management  of  its  fleet  of  F-111  aircraft.  As  a  part  of  the  fatigue  life  enhancement  of 
the  F-111,  the  original  manufacturer  developed  a  testing  procedure  to  detect  cracks  in 
components  of  the  aircraft  that  were  made  from  D6ac  steel.  D6ac  is  a  high  strength 
steel,  and  as  a  consequence  has  an  extremely  low  critical  crack  length.  This 
characteristic  poses  problems  with  regards  to  the  fatigue  life  monitoring  of  the  aircraft. 
Many  critical  areas  of  the  F-111  are  difficult  to  access  with  non-destructive  mspection 
(NDI)  equipment,  making  it  impractical  to  detect  cracking  by  normal  inspection 
methods.  Cold  proof  Load  Testing  (CPLT)  is  the  testing  procedure  developed  to 
overcome  the  complications  presented  by  the  F-111  and  D6ac  steel.  Basically  the  CPLT 
consists  of  first  cooling  the  aircraft  to  -40°F,  where  the  fracture  toughness  of  D6ac  steel 
is  significantly  lower  than  at  room  temperature,  and  then  applying  a  Hirnt  loading 
cycle  to  the  aircraft.  If  the  aircraft  survives  CPLT  it  is  considered  safe  for  a  further 
period  of  time,  before  any  cracks  present  during  CPLT  grow  to  significant  length.  After 
this  period  of  time  the  aircraft  undergoes  another  CPLT  program. 

An  tmdesirable  side-effect  of  the  CPLT  program  is  that  compressive  yielding  at  stress 
concentrators  in  the  upper  plate  of  the  wing  pivot  fitting  (WPF)  leaves  localised  tensile 
stresses  which  facilitate  the  initiation  and  growth  of  fatigue  cracks  under  subsequent 
service  loading.  One  such  area  of  particular  concern  is  that  of  Fuel  Flow  Vent  Hole 
(FFVH)  #13.  The  RAAF  manages  the  structural  integrity  of  its  F-111  fleet  through 
safety  by  inspection,  with  inspection  intervals  derived  from  durability  and  damage 
tolerance  analyses  (DADTA).  The  DADTA  on  which  the  current  inspection  interval  for 
FFVH#13  is  based  took  some  account  of  the  residual  stresses  by  imposing  a  single 
CPLT  load  at  the  start  of  crack  growth  calculations.  The  RAAF  now  requires  a  DADTA 
of  FFVH#13  to  be  performed  incorporating  an  accurate  representation  of  the  residual 
stresses.  AMRL's  role  is  to  determine  the  material  response  for  input  to  this  DADTA, 
which  is  conducted  under  contract  from  the  RAAF.  The  work  described  in  this  report 
forms  part  of  that  broader  AMRL  work  programme  of  which  an  overview  is  given 
in[l]. 

Recent  full  scale  wing  test  data  [2]  indicate  that  the  lower  inboard  comer  of  FFVH#13, 
where  instances  of  cracking  have  occurred,  experiences  strain  levels  ranging  from 
18,000  to  22,000  microstrain  (ps).  The  yield  strain  of  D6ac  is  approximately  ZOOOps,  so 
this  illustrates  the  high  degree  of  plasticity  being  exhibited  by  this  region  of  the 
structure  at  peak  load.  To  accurately  analyse  the  effects  of  such  high  strain  levels  under 
CPLT,  a  Finite  Element  Model  (FEM)  of  the  FFVH#13  region  is  being  employed.  To 
predict  the  residual  stress  left  after  a  CPLT  application,  a  non-linear  FEM  solution 
must  be  performed.  Given  the  cyclic  nature  of  the  CPLT  application,  and  the 
inadequacy  of  classical  plasticity  models  in  predicting  inelastic  cyclic  response,  a 
constitutive  model  was  employed  to  describe  the  inelastic  behaviour  of  the  material. 
This  constitutive  model  applies  a  state  variable,  rate  dependent  approach  which  allows 
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accurate  modelling  of  the  deformation  history  of  the  material  experienced  during 
cyclic  loading.  The  Finite  Element  package  being  used  for  this  analysis  is  PAFEC, 
which  has  the  bruit  in  ability  to  analyse  non-linear  material  problems.  The  standard 
material  models  in  PAFEC  only  make  use  of  classical  techniques  for  modelling 
plasticity,  however  PAFEC  does  allow  the  capability  of  applying  user  written  routines 
such  as  alternative  plasticity  models. 

State  variable  constitutive  modelling  has  been  used  to  improve  predictions  of  the  cyclic 
stress  condition  in  materials.  The  state  variable  approach  to  the  modelling  of  cyclic 
plasticity  attempts  to  model  the  micro  structural  deformation  of  the  material 
experienced  during  the  application  of  cyclic  loading.  This  attribute  results  in  a  more 
accurate  prediction  of  the  material  behaviour.  The  constitutive  model  being  employed 
for  this  work  is  based  on  that  created  by  Ramaswamy,  Stouffer  and  Bodner  [3],  which 
has  been  subsequently  jointly  enhanced  by  AMRL,  Melbourne  University  and  Monash 
University  over  recent  years.  [4,5,6] 

To  develop  the  material  parameters  for  the  initial  prediction  of  the  state  variable 
inelastic  constitutive  model,  sets  of  uniaxial  experimental  data  for  D6ac  steel  were 
required.  This  report  describes  a  series  of  uniaxial  tension/ compression  tests 
performed  to  characterise  the  material  behaviour  of  D6ac  steel  under  cyclic  inelastic 
loading.  These  data  sets  were  then  used  to  determine  the  initial  values  of  the  unified 
constitutive  model's  material  constants.  This  first  estimate  was  refined  using  a 
FORTRAN  computer  program  which  represents  the  uniaxial  testing  condition  using 
the  unified  constitutive  equations. 

In  a  parallel  program  of  work,  the  theoretical  and  analytical  solution  techniques  of  the 
unified  constitutive  model  were  developed  in  references  3,  5  and  6.  The 
implementation  of  the  analysis  software  into  the  PAFEC  finite  element  program  can  be 
fotmd  m  references  4  and  7.  The  achievements  of  this  work  now  make  it  possible  to 
perform  the  required  plasticity  analysis,  accurately  and  efficiently.  This  capability  is 
being  used  on  the  F-lllC  FFVH#13  problem  to  solve  for  the  large  cyclic  plastic  strains 
that  occur  during  the  CPLT. 


2.  Theory 


The  unified  constitutive  model  developed  by  Ramaswamy,  Stouffer  and  Bodner  [3] 
consists  of  the  inelastic  flow  equation  and  the  rate  forms  of  the  state  variables;  back 
stress  Qi|,  and  drag  stress  Z.  The  drag  stress,  Z,  is  taken  as  a  scalar  value  as  its 
evolution  (interaction  of  dislocations  with  subgrains,  cells,  precipitates,  etc)  is  largely 
orientation  independent  [8].  Equations  for  the  multiaxial  form  of  the  tinified 
constitutive  model  are: 


8^.  =  Dexp 


(V 


2 


DSTO-TR-0556 


where  sj-  is  the  inelastic  strain  rate,  K2  is  the  second  invariant  of  the  inelastic  strain 

rate  tensor,  Sij  is  the  deviatoric  stress  tensor,  and  D  and  n  are  temperature  and  strain 
rate  dependent  material  parameters.  For  D6ac  steel,  at  room  temperature,  D=10000  for 
tensile  and  fatigue  responses  (for  strain  rates  up  to  1.0  sec-^),  and  n=3  for  strain  rate 
insensitivity.  Here  the  evolution  equation  for  the  back  stress  Qjj  is  defined  as: 

where  f2  is  a  material  parameter  and  Sij  is  the  deviatoric  stress  tensor.  The  inelastic  back 
stress  Q  |j  is  calculated  by  integrating  the  following. 


Q 

ij 

Q 


'eff 


(3) 


where  Q  |j  is  the  inelastic  back  stress  rate,  Qmax  is  the  maximum  value  of  back  stress 
expected,  is  the  effective  inelastic  strain  rate  and  fi  is  a  hardness  related  material 
parameter.  The  effective  inelastic  strain  Sgg,  is  given  by. 


where  s|j  is  the  inelastic  strain  tensor. 


The  drag  stress  Z  is  given  by, 

Z  =  +  (Zq  -  Z-i)  exp  ( -  ms'g)  (5) 

where  Zo  is  the  initial  value  of  drag  stress,  Zi  is  the  final  value  of  drag  stress  and  m 
defines  the  rate  at  which  Z  evolves  from  Zo  to  Zi. 

From  the  above  set  of  equations,  the  following  characteristic  constants  were  required 
to  be  derived  from  uniaxial  test  data:  Do,  n,  fp  f2/  ^max'  ^0'  ^1'  determination 

of  these  parameters  are  dependent  on  a  number  of  mathematical  relationships  which 
describe  the  state  variables  in  terms  of  material  characteristics  seen  in  typical  xmiaxial 
tests.  These  relationships  are  related  to  the  following  material  characteristics: 
saturation  of  a  tensile  response  curve,  the  Bauschinger  effect,  stress  relaxation  and 
strain  recovery. 

The  uniaxial  form  of  equations  1  to  4  are  utiHsed  in  section  7.1  to  provide  initial 
estimates  of  the  constitutive  state  variables,  which  are  then  improved  upon  with  the 
use  of  a  perturbation  computer  program  to  arrive  at  the  final  set  of  values,  described  in 
section  9. 
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3.  Test  Specimens 

The  specimen  design  employed  for  this  test  series  was  modified  from  that  used 
previously  at  AMRL  [9].  These  specimens  incorporated  threads  at  both  ends  so  that 
they  could  be  screwed  into  a  set  of  custom  made  holding  blocks  and  then  inserted  into 
standard  testing  machine  grips,  (see  Figure  1)  The  specimens  were  cylindrical  in 
section  and  127nun  in  length  (including  threaded  ends)  with  a  9.5mm  diameter  test 
section  (see  Figure  2).  There  were  10  specimens  in  total,  with  each  specimen  having 
been  previously  heat  treated  to  a  level  consistent  with  that  of  the  WPF. 

It  was  noted  that  some  of  these  specimens  were  slightly  distorted  as  a  result  of  the 
manufacturing  and  heat  treatment  processes.  This  subsequently  affected  the  alignment 
process  (see  Section  8.1). 

In  order  to  obtain  the  most  accurate  set  of  results,  the  alignment  of  the  specimens  in 
the  testing  machine  was  paramount.  Each  specimen  was  fitted  with  a  set  of  four  strain 
gauges  at  ninety  degree  intervals  aroxmd  the  test  section.  The  strain  gauges  were 
employed  to  aid  the  alignment  of  the  specimens  whilst  being  gripped  in  the  testing 
machine.  They  were  applied  using  high  temperature  cure  adhesives  to  provide  the 
highest  bond  failure  strain.  The  strain  gauges  were  expected  to  drift  xmder  the 
application  of  high  strains  and  their  output  would  be  unreliable  once  the  strains 
reached  these  high  levels  (greater  than  10,000  ps) [10]. 

4.  Test  Equipment 

The  following  testing  equipment  was  used  throughout  the  test  series. 

-  Instron  250  kN  universal  servo  hydraulic  testing  machine  located  in  the  AMRL 

Fatigue  Laboratory. 

-  Instron  Extensometer  (Cat  no.  2620-602),  25mm  gauge  length, 

±10%  strain  range. 

-  Instron  Extensometer  (Cat  no.  2620-603),  10mm  gauge  length, 

±10%  strain  range. 

-  Mitutoyo  Digmatic  Indicator  (Dial  Gauge)  (Code  no.  543-612), 

Accuracy  0.02iiun  @  20°C. 

The  25mm  extensometer  was  placed  along  the  axis  of  the  specimen  (see  Figure  1)  for 
all  of  the  tests.  The  10mm  extensometer  was  also  positioned  across  the  centre  of  the  test 
section,  however  it  was  placed  135  degrees  around  the  circrunference  of  the  test  section 
with  respect  to  the  25mm  extensometer  for  specimens  3, 6  and  7  only. 

Test  data  were  recorded  on  a  PC,  using  the  in-house  program  ISGAR  version  1.85,  with 
each  of  the  following  parameters  measured: 
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-  Time 

-  Load  (raw  voltage) 

-  Extensometer  reading  (raw  voltage) 

-  Four  strain  gauge  readings  (raw  voltage) 

-  Testing  machine  LVDT  (for  initial  tests  only) 


Also  during  the  testing  the  following  real  time  plots  were  recorded: 

-  Load  and  Extensometer  strain  vs  time  on  chart  recorder 

-  Load  vs  Extensometer  strain  on  HP  XY  plotter. 


5.  Test  Procedure 


5.1  Alignment  Procedure 

Prior  to  commencement  of  the  testing  program  each  specimen  was  inspected  for  any 
possible  distortion  due  to  the  heat  treatment  process.  Most  of  the  specunens  were 
foxmd  to  be  slightly  bowed,  with  some  showing  a  more  marked  effect  than  others. 

Installation  of  the  specimens  into  the  testing  machine  began  with  the  specimen  being 
screwed  into  the  custom  made  holders  at  each  end,  and  the  locking  nuts  attached,  but 
not  tightened  (see  Figure  1).  The  next  step  was  to  clamp  the  specimen  into  the  top  grip 
of  the  machine,  leaving  the  bottom  end  free.  By  applying  a  dial  gauge  to  the  edge  of 
the  top  holder  the  specimen  could  be  approximately  aligned  in  the  vertical  direction  by 
continually  readjusting  the  top  holder.  At  this  point  the  bottom  end  of  the  specimen 
was  clamped  and  the  nuts  were  tightened.  The  strain  gauge  values  were  recorded  and 
a  small  compressive  strain  of  -2000  ps,  was  applied.  Any  bending  of  the  specimen  was 
observed  through  the  tmequal  response  of  the  four  gauges. 

Next,  both  nuts  were  loosened  and  the  lower  specimen  holder  was  then  shimmed  in 
the  appropriate  direction  to  remove  any  apparent  bending  observed.  The  process  was 
then  repeated  until  any  observed  bending  effects  were  minimised.  This  method  was 
time  consuming  but  accurate,  considering  that  the  specimens  already  possessed  a 
certain  degree  of  bend.  This  method  of  aligning  the  specimens  also  yielded  an  accurate 
estimate  of  the  degree  of  bending  present  within  a  specimen,  thus  providing  a  greater 
xmderstanding  of  &e  results. 

The  final  two  tests,  specimens  6  and  7,  were  aligned  in  a  different  manner  due  to 
experience  gained  from  the  previous  tests.  These  specimens  were  first  inserted  into  the 
machine  and  the  nuts  tightened.  A  small  compressive  strain  was  applied  and  the  strain 
gauge  results  noted.  The  specimen  was  then  removed  from  the  machine  and  then 
rotated  90  degrees  within  the  holders  and  then  replaced.  Again  a  small  compressive 
strain  was  applied  and  the  results  noted.  This  was  repeated  a  final  time  at  a  rotated 
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angle  of  45  degrees  from  the  initial  position.  From  these  three  sets  of  results  the 
bending  axis  was  calculated.  The  hydrauUc  grips  used  only  allow  alignment  in  one 
plane.  The  specimen  was  then  installed  so  that  its  bending  axis  was  in  the  plane  of  the 
testing  machine.  Specimen  distortion  is  then  minimised  by  shimming  the  holder  in  the 
hydraulic  grips. 


5.2  Specimen  Loading  Configuration 

The  strain  ranges  for  the  tests  were  chosen  to  simulate  the  strain  levels  seen  in  the 
FFVH#13  region  at  100%  CPLT  load  as  given  in  reference  [9].  A  symmetric  tension  / 
compression  load  cycle  was  used  to  characterise  the  material,  whereas  a  non- 
symmetric  CPLT  cycle  was  used  to  represent  the  actual  strain  cycle  experienced  by  the 
FFVH#13  region. 

Some  materials  have  shown  a  non-symmetrical  tension  and  compression  monotonic 
loading  resporise  [11].  This  phenomenon  was  investigated  in  this  series  of  tests.  Most 
uniaxial  data  is  tension  oriented  only.  However,  the  CPLT  cycle  involves  the  initial 
application  of  large  compressive  loads.  Therefore  it  was  necessary  to  investigate  both 
the  tension  and  compression  loading  regimes. 

Previous  work  has  shown  that  D6ac  steel  is  strain  rate  independent.  Incorporated  into 
the  test  program  was  the  validation  of  this  material  characteristic.  The  strain  rate  of 
60ps/sec  is  representative  of  the  actual  CPLT  application  strain  rate  and  therefore  this 
was  used  as  the  baseline. 

To  fully  characterise  the  material,  the  unified  constitutive  model  requires  a  knowledge 
of  stress  and  strain  relaxation.  In  this  test  series  only  the  stress  relaxation  was 
investigated,  in  order  to  evaluate  the  back  stress  state  variable,  Q  (see  Section  7.2). 

The  following  is  a  detailed  description  of  the  applied  loading  and  a  description  of  each 
test  specimen. 

5.2.1  Specimen  No  1:  Basic  tension/ compression  cycle  at  60ps/ sec 

Under  extensometer  strain  control,  cycle  to  ±20,000  ps  tension/ compression  at  a  strain 
rate  of  60  ps/sec,  beginning  with  the  tension  cycle  (see  Figure  3).  This  test  was  to 
investigate  the  symmetry  of  a  single  inelastic  strain  loop. 

5.2.2  Specimen  No  2:  Basic  compression/ tension  cycle  at  60ps/ sec 

Repeat  test  of  specimen  nxunber  1,  however  the  cycle  began  with  a  compression  cycle. 
Again  the  test  was  imder  strain  control  and  at  60  ps/ sec  (see  Figure  4).  The  main  aim 
of  this  test  was  to  compare  the  loop  shape  of  a  single  cycle  beginning  with  a 
compression  application  as  opposed  to  that  beginning  with  a  tension  load. 
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5.2.3  Specimen  No  3:  Investigation  of  specimen  bending  effects 

This  test  was  used  to  monitor  the  possible  bending  effects  seen  during  the  previous 
test  sequences.  Initially  the  specimen  underwent  a  number  of  separate  20,000ps/- 
20  OOOps  cycles  with  two  extensometers  of  differing  gauge  lengths  applied  to  the  test 
section.  After  these  cycles,  a  number  of  tests  at  increasing  symmetric  strain  levek  were 
conducted  to  quantify  the  monotonic  saturated  stress  state  for  D6ac  steel,  i.e.  the  load 
and  strain  levels  when  the  stress/ strain  curve  becomes  flat.  Strain  cycles  were  applied 
up  to  35,000|J.s/-35,000|.is,  in  increments  of  ±5,000|is  xmtil  the  specimen  failed  m 
compression.  This  failure  was  a  direct  result  of  the  specimen  bending  out  of  plane  on 
the  compression  cycle. 

5.2.4  Specimen  No  4:  Basic  tension/ compression  cycle  at  lOpe/sec 

Under  strain  control,  cycle  to  ±20,000  pe  tension/ compression  at  a  strain  rate  of  10 
ps/sec,  beginning  with  the  tension  cycle  (see  Figure  5).  This  test  was  used  to 
investigate  the  strain  rate  sensitivity  of  D6ac  steel. 

5.2.5  Specimen  No  5:  Basic  tension/ compression  cycle  at  600pe/ sec 

Under  strain  control,  cycle  to  ±20,000  pe  tension/ compression  at  a  strain  rate  of  600 
ps/sec,  beginning  with  the  tension  cycle  (see  Figure  6).  This  test  was  used  to 
investigate  the  strain  rate  sensitivity  of  D6ac  steel. 


5.2.6  Specimen  No  6;  Stress  relaxation  under  strain  control 

This  test  began  by  applying  a  tensile  loading  imder  strain  control  at  the  rate  of  60 
ps/sec.  At  15,000  the  strain  was  held  for  2  horns  and  data  recorded  initially  every 
second  for  3  minutes,  then  at  20  second  intervals.  From  this  point  the  specimen  was 
loaded  up  to  20,000  ps  and  then  down  to  -15,000  Here  the  strain  was  held  for  2  hours 
and  the  stress  relaxation  was  recorded  for  the  compression  case  every  second  for  3 
minutes  then  at  20  second  intervals.  The  strain  was  continued  to  -20,000  ps  and  then 
the  specimen  was  cycled  between  20,000  ps  and  -20,000  ps,  until  cyclic  saturation,  (ie. 
no  change  in  stress/ strain  loop  observed  between  subsequent  cycles)  This  took 
approximately  five  complete  cycles.  The  above  process  was  repeated  three  times, 
taking  approximately  5  cycles  for  each  loop  saturation. 

This  test  was  specific  to  the  calculation  of  the  constitutive  equation  state  variables. 
Although  die  constitutive  relationships  realistically  only  require  stress  relaxation  data 
on  the  first  cycle  for  their  development,  this  test  provided  data  on  any  differences  in 
the  stress  relaxation  levels  after  applying  a  cyclic  loading. 

5.2.7  Specimen  No  7;  Stress  relaxation  under  strain  control 
Repeat  test  of  specimen  number  6. 


7 


DSTO-TR-0556 


5.2.8  Specimen  No  8: 100%  CPLT  representative  cycle 

Apply  CPLT  cycle  (2  cycles  of  Figure  7)  at  a  constant  strain  rate  of  60  ps/ sec. 

5.2.9  Specimen  No  9:  75%  CPLT  representative  cycle 

Apply  75%  CPLT  cycle  (2  cycles  of  Figure  8)  at  a  constant  strain  rate  of  60  ps/  sec. 

5.2.10  Specimen  No  10;  50%  CPLT  representative  cycle 

Apply  50%  CPLT  cycle  (2  cycles  of  Figtue  9)  at  a  constant  strain  rate  of  60  ps/sec. 


6.  Results  and  Discussion 


The  results  of  this  experimental  test  program  were  used  to  determine  the  state 
variables  for  the  unified  constitutive  model  (see  Section  7.1).  The  raw  data  from  each 
specimen  are  presented  in  the  form  of  load  versus  microstrain  curves  (Figures  10-19), 
the  microstrain  outputs  are  from  the  25mm  extensometer  unless  otherwise 

noted. 

Observations  made  during  the  testing  program  are  discussed  in  detail  in  the  following 
sections; 


6.1  Specimen  Alignment 

As  mentioned  previously,  the  specimens  were  aU  slightly  distorted,  which  made 
accurate  specimen  alignment  extremely  difficult  from  the  outset.  The  specimens  were 
aligned  as  accurately  as  possible  to  ensure  bending  was  kept  to  a  minimum.  Most  of 
the  tests  involved  only  one  or  two  cycles  and  therefore  it  could  be  assumed  that 
bending  effects  never  progressed  to  a  level  which  was  detrimental  to  die  results. 
Specimen  3  however,  which  involved  multiple  cycle  testing,  proved  susceptible  to 
bending.  This  can  be  illustrated  by  Figure  12,  where  specimen  3  experienced  bending 
effects  after  it  had  imdergone  a  number  of  continuous  cycles.  Specimen  3  made  use  of 
two  extensometers  of  differing  gauge  lengths  placed  at  different  circumferential 
positions  arotmd  the  working  section.  One  extensometer  of  larger  gauge  length 
(25mm)  controlled  the  experiment,  whilst  the  second  extensometer  of  smaller  gauge 
length  (10mm)  was  also  placed  across  the  centre  of  the  specimen  test  section,  but  at  135 
degrees  aroimd  from  the  larger  one.  No  significant  bending  occurred  within  the  first 
few  cycles,  however  after  10  cycles  bending  of  the  specimen  became  obvious.  This 
effect  grew  as  the  ntunber  of  cycles  increased.  Specimens  6  and  7,  which  also  involved 
multiple  cycle  tests,  did  not  produce  the  same  degree  of  bending  as  did  specimen  3. 
The  bending  observed  for  these  specimens  was  extremely  minimal,  as  both  attached 
extensometers  returned  the  same  strain  outputs. 
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OveraU,  this  confirmed  that  the  degree  of  bending  was  specimen  dependent,  and  that 
specimen  3  was  the  only  specimen  to  exhibit  the  effects  of  bending  to  a  significant 

degree. 


6.2  Uniaxial  Specimen  Design 

From  the  beginning  of  the  testing  program  it  was  recognised  that  the  specimen  design 
was  not  optimum  for  high  tension  and  compressive  testing.  The  ideal  aspect  ratio  of 
these  types  of  specimens  should  be  approximately  1.5  [12],  which  would  result  m 
minimal  bending  and  bulging.  Bulging  was  not  observed  in  any  of  the  specimens 

tested. 

The  application  of  the  load  to  the  test  section  could  also  be  improved.  The  current 
setup  entails  the  specimen  being  screwed  into  a  set  of  high  strength  tool  steel  holdmg 
blocks  and  then  inserted  into  the  test  machine  grips.  The  load  is  transferred,  m  shear, 
through  the  threaded  section  of  the  specimen.  Before  the  testing  commenced  it  was 
thought  that  the  threads  on  the  specimens  might  deform  under  loadmg  and 
consequently  result  in  an  inefficient  appHcation  of  load  to  the  test  section.  As  the  test 
program  progressed  it  was  found  that  the  specimen  holders  threads  were  deforming. 
This  was  discovered  through  the  increasing  difficulty  in  inserting  the  specimens  mto 
the  holders  as  the  tests  progressed.  It  is  thought  that  this  deformation  could  be  a  result 
of  the  bending  of  Ihe  specimens  implying  uneven  loads  through  the  threads  of  the 
specimen  holders,  and  thus  resulting  in  thread  deformation.  Considering  the  extent  of 
thread  deformation  experienced,  it  would  be  beneficial  if  an  alternative  gripping 
technique  could  be  devised  for  future  tests. 

In  order  to  determine  that  the  test  section  would  exhibit  a  uniform  stress  regime 
throughout  the  cross  section,  an  idealised  elastic  FEM  analysis  of  Ihe  specimen  was 
performed.  This  analysis  is  detailed  in  Appendix  A  and  shows  that  the  specimens 
experienced  a  constant  stress  application  across  the  working  section. 


6.3  Cyclic  Saturation 

An  interesting  material  characteristic  seen  throughout  all  of  the  multiple  cycle  tests 
was  the  way  in  which  the  material  cyclically  saturated.  Most  materials  take  a  few 
cycles  to  cyclically  saturate,  however  D6ac  steel  is  cyclicaHy  saturated  after  the  first 
cycle.  This  can  be  seen  in  Figures  15  and  16  where,  once  the  material  initiaUy  yields  on 
the  first  cycle,  the  subsequent  cycles  lie  over  one  another.  This  effect  is  consistent 
throughout  all  of  the  multiple  cycle  tests,  tension,  compression  and  CPLT. 
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6.4  Inelastic  Saturation 

Another  interesting  material  characteristic  observed  during  the  testing  was  the  high 
level  of  strain  required  to  reach  inelastic  saturation  (when  the  stress/ strain  curve 
becomes  flat).  All  tests  to  20,000  |xs  showed  that  the  material  had  not  inelastically 
saturated.  Furthermore,  specimen  3  which  reached  approximately  35,000  ps  before 
failxu-e,  clearly  did  not  demonstrate  that  the  material  had  inelastically  saturated.  From 
these  data  it  could  be  estimated  that  D6ac  steel  would  inelastically  saturate  at  around 
40,000  ps. 


6.5  Stress  Relaxation 

In  addition  to  the  two  separate  tests  specifically  conducted  to  investigate  the  stress 
relaxation  properties  of  D6ac  steel,  stress  relaxation  could  be  seen  during  every  test. 
Strain  control  was  used  for  all  the  tests,  and  as  a  result  of  this,  stress  (or  load) 
relaxation  was  observed  when  a  specimen  was  returned  to  zero  strain  at  the  end  of 
each  test.  All  test  results  showed  the  load  to  relax  approximately  30.27  MPa 
instantaneously  upon  test  completion.  If  the  specimen  was  left  for  about  5  minutes, 
this  relaxation  progressed  to  arotmd  56.43  MPa.  It  is  interesting  to  note  that  the  2  hour 
strain  holds  for  the  stress  relaxation  tests  showed,  on  average,  a  88.88  MPa  load 
relaxation  (see  Figures  15  &  16).  Figure  20  shows  a  section  of  test  data  for  specimen  6 
displayed  in  load  time  format.  The  bulk  of  the  stress  relaxation  occurs  within  5 
minutes  which  is  seen  in  all  of  the  tests. 


6.6  Bauschinger  Effect 

The  Bauschinger  effect,  which  is  the  lowering  of  tiie  absolute  value  of  the  elastic  limit 
in  compression  following  a  previous  tensile  loading  [11],  was  prominent  throughout 
all  of  the  tests,  as  it  was  for  previous  D6ac  steel  tests.  Through  analysis  of  a  basic 
symmetric  tension/ compression  strain  cycle  (Figure  10),  it  can  be  seen  that  D6ac  yields 
quite  sharply  on  the  initial  yield  cycle  and  then  more  gradually  on  subsequent  cycles. 
This  effect  within  itself  leads  to  the  occurrence  of  the  Bauschinger  effect,  as  the 
compression  yield  following  the  initial  tensile  yield  is  vastly  different  from  that 
experienced  within  the  initial  yield.  This  phenomenon  persists  throughout  multiple 
cycle  tests  (Figtues  12,15,16),  as  the  cyclic  loops  are  virtually  identical.  Previous  testing 
also  encotmtered  the  Bauschinger  effect,  which  further  vahdates  its  authenticity  as  a 
characteristic  of  D6ac  steel. 


6.7  Material  Anomalies 

Another  interesting  characteristic  displayed  by  D6ac  steel  during  the  cyclic  testing  was 
the  differing  values  of  Young's  Modulus  (E)  for  initial  loading  and  subsequent  xmload. 
For  a  single  20,000/ -20,000  ps  loop  beginning  with  a  tension  load  (see  Figtire  10),  the 
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Voting's  Modulus  was  the  same  for  the  initial  tension  upload  and  the  unload  after  a 
compressive  load.  However,  a  lower  value  of  E  was  observed  when  unloading  from 
any  tensile  load  (see  Figures  17,18,19).  This  phenomenon  was  quantitatively  consistent 
throughout  aH  of  the  tests  and  was  exhibited  by  both  extensometer  and  stram  gauges 
alike.  Unloading  in  the  elastic  region  of  the  tensile  curve  did  not  produce  a  different 
value  of  Young's  Modulus,  only  in  unloading  after  plastic  deformation  was  this 
characteristic  evident. 

Another  material  anomaly  observed  was  the  differing  maximum  load  values  reached 
in  tension  and  compression  for  a  symmetric  strain  cycle  appHcation.  The  maximum 
load  value  reached  in  compression  was  consistently  10  kN  higher  than  that  seen  in 
tension  for  all  of  the  symmetric  tests.  Given  the  associated  consistency  of  the  differing 
values  of  E,  these  two  characteristics  are  possibly  inter-related. 

There  exist  two  possible  explanations  of  these  material  anomahes.  Firstly,  if  the 
deformation  mechanism  of  D6ac  is  slip  (most  likely  mode  of  deformation)  then  the 
cyclic  cxirve  and  Yotmg's  Modulus  should  be  symmetric  in  tension  and  compression. 
In  this  case  the  specimens  may  be  bending  due  to  an  imperfect  alignment,  thus  causing 
the  differing  value  of  E  and  maximum  load,  however  the  consistency  of  the  results 
across  aU  of  the  specimens  discounts  this  theory.  It  may  be  possible  that  sHp  is  not  the 
only  deformation  mechanism.  Grey  cast  iron  is  not  symmetric  in  tension  ^d 
compression  and  its  deformation  occurs  through  mechanisms  other  than  slip  [13].  It 
could  be  possible  that  D6ac  experiences  mechanisms  other  than  slip  which  results  in  a 
non-symmetric  response.  It  is  also  possible  that  micro  voids  could  have  been 
introduced  into  the  specimen's  structure  during  heat  treatment.  As  the  specimens 
imdergo  plastic  deformation  in  the  tensile  region,  the  micro  voids  would  be  stretched 
open.  Upon  load  reversal  the  specimens  would  effectively  possesses  a  smaller  cross 
sectional  area,  thus  resulting  in  a  decreased  value  of  Young's  Modulus.  As  the 
specimeris  are  compressed,  the  micro  voids  would  be  closed,  resulting  in  the  original 
Young's  Modulus  occurring  when  the  specimen  xinloads  from  the  maximum  tension 
load.  The  modified  Yotmg's  Modulus  does  not  occur  in  the  initial  tension  loading,  as 
the  micro  voids  would  only  expand  out  to  a  significant  amount  during  tensile  plastic 
deformation.  This  explanation  would  require  a  microscopic  analysis  of  the  material  to 
confirm  the  existence  of  these  voids,  and  thus  their  effect  on  the  material's  behaviour. 


7.  Constitutive  Model  State  Variable  Determination 

7.1  Initial  Estimation  of  State  Variables 

It  should  be  noted  that  the  equations  given  in  this  section  are  used  as  simplified 
asymptotic  approximations  only.  Their  purpose  is  to  provide  a  first  estimate  of  the 
values  of  the  state  variables,  so  that  they  can  be  refined  to  yield  an  accurate  prediction 
of  the  material  response.  Once  such  an  initial  estimate  of  the  state  variables  was 
obtained  they  were  input  into  a  FORTRAN  computer  program  which  contained  the 
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unified  constitutive  equations  1  to  5.  The  program  was  designed  to  represent  the 
uniaxial  test  condition  under  strain  control  and  it  allowed  the  perturbation  of  the  state 
variables  and  the  subsequent  analysis  of  the  response  compared  to  uniaxial  test  data. 


7.1.1  Tensile  Response  Data 

The  values  of  Do  and  n  have  been  selected  based  upon  previous  experience  for  strain 
rate  independent  materials  [3,8].  At  tensile  saturation,  the  back  stress  is  at  a  maximum, 
termed  Qmax,  whilst  the  drag  stress  is  near  it's  initial  value,  Zo,  and  remains  almost 
constant.  For  isothermal  conditions  the  temperature  can  be  omitted  from  equation  1. 
Therefore  for  a  saturated  uniaxial  tensile  loading  (where  a  is  at  saturation  and  is  a 
maximxim),  equation  1  can  be  re-written  in  the  following  form; 


=  ^P„exp 


■  If 
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where  the  term,  cr  -  Q/|<t  -  Q|  in  equation  1  is  equal  to  +1,  and  osat  is  the  value  of  the 
saturated  stress.  Inverting  equation  6  yields  the  following  which  can  be  used  for 
obtaining  a  relationship  between  flmax  and  Zo. 
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Here  el  is  the  initial  inelastic  strain  rate.  For  a  constant  strain  rate  application  Sq  is 
equal  to  the  overall  strain  rate.  Therefore  from  specimen  1  (Figure  10),  Sp  =  0.00006 
s/sec.  Substituting  this  and  the  relevant  values  of  Do  and  n  into  equation  7  yields: 
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=  1.83 


which  leads  to  the  following  equation: 

Z.=1.83(<T„,-n^,)  (8) 

Equation  8  is  used  later  in  conjtmction  with  equations  11,12  and  16  to  solve  for  Qmax 
and  Zo. 
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7.1.2  Bauschinger  Effect 

Through  analysis  of  the  flow  equation  1,  it  can  be  seen  that  the  over  stress  (resultant 
stress,  a-Q)  to  produce  yield  in  tension  must  be  the  same  as  the  over  stress  to  produce 
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yield  in  compression,  for  constant  values  of  n.  Do  and  Zo.  Therefore  for  the  first  cycle  of 
a  fatigue  test  when  Z  =  Zo,  equation  1  can  be  rewritten  in  the  following  form: 

(9) 

where  Oyt  and  Gyc  are  the  values  of  yield  stress  for  tension  and  compression 
respectively.  Through  analysis  of  Figure  21  the  following  two  relationships  between  £2 
and  Qm  can  be  discerned,  (see  [8]) 


(10) 

^yc  =  ^max  -  f 2' 

(11) 

From  specimen  1  data  the  following  results  were  extracted: 

0yt  =  1340  MPa 
Oyc  =  -141  MPa 
Aa  =  1763  MPa 

Inserting  these  values  into  equation  11  leads  to  the  following  relationship: 

-141xl0‘=a^-/4l763xl0‘)  (12) 


7.1.3  Saturation  Response 

Next  we  make  use  of  the  parameter  fi  which  describes  the  inelastic  strain  to  reach 
saturation,  see  [8]. 
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where  is  the  accumulated  inelastic  strain  at  the  onset  of  saturation. 

7.1.4  Stress  Relaxation 

Stress  relaxation  properties  can  be  made  use  of  in  the  following  fashion  (see  [8])  to 
relate  ii  to  Dmax  (see  Figure  22). 


f2  = 


c, 


^max 


\AaAs\  +  C,a^ 


(14) 


where 
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+  (15) 

Es  sat 

From  specimeir  6  the  foUowing  data  were  determined; 

GB  =  1445  MPa  OA  =  1535  MPa 
ayt  =  1340  MPa  E  =  203  GPa 

4,  =40,000 

Subsatutmg  these  values  into  equations  14  and  15  yields  the  following  equation: 

(1.0252)n„  -  (1429  X  10‘)/,  =  107  x  10‘  (16) 

Therefore  solving  equations  12  and  16  for  ii  and  Qmax: 

amax=  1031x10"  MPa 
f2  =  0.665 

Substituting  Qmax  into  equation  13  yields; 
fi  =  118.565  GPa 

And  finally  substituting  Qmax  into  equation  8; 

Zo=  1081.53x10"  MPa 

These  values  yield  an  initial  estimate  of  the  material  constants  for  the  constitutive 
model. 


7.2  Refinement  of  Initial  Predictions 


The  initial  values  of  the  material  constants  for  D6ac  steel  were  refined  by  coinpanng 
predictions  with  symmetric  cyclic  strain  data  as  weU  as  CPLT  cycle  data. 
This  produced  the  final  set  of  constitutive  constants  listed  below. 


=  1250  MPa 
Zo  =  1081  MPa 

fq  =  135  GPa 

f2  =  0.585 


The  above  set  of  values,  are  relatively  different  to  the  initial  estimates  calcukted  using 
ctl^totiva  relatiorwhips.  The  vaL  of  D™.  is  larger  by  about  200  MPa.  This  rs  due 
to  the  different  maximum  loads  seen  during  a  smgle  tension/ compression  cycle,  s 
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mentioned  in  the  uniaxial  test  data  discussion,  the  maximum  tension  load  was 
consistently  about  10  kN  lower  in  comparison  to  the  maximum  compression  load.  The 
model  does  not  allow  for  different  maximum  tension  and  compression  loads  and 
therefore  a  higher  value  of  Qmax  was  required  to  be  able  to  predict  an  average 
tension/ compression  response. 

Also  the  values  of  fi  and  fz  are  slightly  modified  with  respect  to  their  initial 
estimations.  Basically,  these  variables  control  the  point  at  which  the  material  yields 
and  also  the  approximate  strain  at  which  the  material  saturates.  By  analysing  a 
symmetric  strain  cycle  for  D6ac  steel  (Figure  10)  it  can  be  seen  that  the  material  has  a 
very  sharp  initial  yield  response.  However,  after  the  material  has  initiaUy  yielded,  the 
subsequent  cyclic  response  is  quite  smooth.  Therefore,  to  be  able  to  predict  this  vast 
difference  in  the  initial  yield  and  subsequent  yields,  the  values  of  fi  and  ii  were 
modified  accordingly.  The  variable  m,  in  effect  controls  the  cycHc  softenmg  of  the 
material.  The  final  value  of  120  was  obtained  by  varying  m  within  the  perturbation 
program  until  an  acceptable  response  was  achieved. 

Results  of  the  above  state  variables  can  be  seen  compared  to  experimental  data  in 
Figures  23  and  24. 


7.3  Discussion  of  Constitutive  Model  Prediction 

As  can  be  seen  from  Figures  23  and  24  the  numerical  prediction  of  the  material 
behaviour  of  D6ac  steel  is  relatively  good,  however  the  prediction  is  erroneous  around 
certain  areas  of  the  curve.  Examining  Figure  23,  which  depicts  a  symmetric  20,000  ps/- 
20,000  ps  loop,  the  largest  discrepancy  of  the  numerical  prediction  is  seen  to  occxir 
within  the  reverse  yield  area,  with  a  maximum  difference  occurring  at  around  -5000  ps. 
This  large  difference  can  be  attributed  to  the  reverse  yielding  predicted  by  the 
numerical  model  after  the  compression  cycle,  compared  with  the  corresponding 
reverse  yield  Yoxmg's  Modulus  displayed  by  the  experimental  results.  As  explained 
earlier  the  gradient  of  the  elastic  unload  from  tension  yield  is  different  from  that 
exhibited  within  initial  loading  and  unloading  from  compression  yield.  This  is  an 
uncommon  characteristic,  and  the  numerical  model  is  incapable  of  predicting  the 
correct  response.  It  is  also  apparent  from  Figure  23  that  the  model  has  difficulty 
predicting  the  strain  hardening  of  the  material  after  the  initial  yield.  The  model 
predicts  a  more  gradual  yielding  than  that  seen  in  the  uniaxial  tests.  These  limitations 
coupled  with  the  differing  values  of  Pmax  and  Pmin,  make  it  impossible  to  obtain  an 
exact  prediction  of  the  material  response. 

The  constitutive  model  is  based  theoretically  upon  the  development  of  the  material 
microstructure  consistent  with  the  physical  properties  of  the  material.  The  current 
formulation  of  the  model  is  limited  in  its  capability  to  model  the  observed 
inconsistencies  in  the  material  response.  This  limitation  would  occur  with  other 
plasticity  models,  as  the  material  anomalies  displayed  by  the  uniaxial  tests  do  not 
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agreG  with  the  expected  material  behaviour,  and  have  not  been  considered  in 
formulating  the  models. 

Another  interesting  physical  characteristic,  which  leads  to  the  inaccurate  prediction  of 
the  material  response,  is  the  sharp  yielding  exhibited  by  D6ac  steel  on  the  first  yield,  as 
compared  with  subsequent  yields,  as  can  be  seen  from  Figure  23.  This  is  hard  for  the 
constitutive  model  to  predict,  because  the  constitutive  model  is  based  on  a  single  set  of 
state  variables  which  remain  constant  throughout  the  evolution  of  the  material 
behaviour.  This  excludes  the  variables  Zi  and  Zo  which  describe  the  upper  and  lower 
boimds  of  the  drag  stress,  and  evolve  throughout  the  material  deformation.  Other 
variables,  such  as  fi  and  fi,  for  example,  which  describe  the  evolution  of  the  back  stress, 
Qm,  do  not  evolve. 

Basically,  the  terms  fi  and  fa  control  the  yield  point  of  the  material.  In  general,  the 
larger  the  value  of  fi  the  higher  the  stress  reached  before  the  onset  of  yield,  whilst  the 
higher  the  value  of  fa,  the  lesser  the  value  of  strain  whereby  the  material  saturates.  The 
opposite  effects  are  also  true.  If  fi  and  fa  were  the  only  controlling  mechamsms  of  the 
material  yield,  then  the  same  shape  of  yield  for  the  initial  loading,  and  for  subsequent 
cyclic  loading  would  be  experienced,  however  this  is  not  the  case.  The  hardening 
parameters  fi  and  fa  effectively  combine  with  the  evolving  value  of  the  drag  stress,  Z, 
which  thus  leads  to  the  alleviation  of  the  sharp  yield  effect  seen  for  the  initial  loading. 

However,  there  seems  to  exist  a  Umit  to  which  these  parameters  can  combine  to 
produce  a  smooth  transition  from  uutial  yield  to  cyclic  yield.  As  can  be  seen  in  Figures 
23  and  24,  the  constitutive  model  cannot  predict  the  initial  yield  point  and  then 
accurately  predict  the  extremely  gradual  yield  depicted  by  the  experimental  data  for 
the  cycHc  loading.  This  is  because  the  only  variable  that  evolves  from  the  initial  yield 
to  the  cyclic  yield  is  the  drag  stress.  The  values  of  fi  and  remain  constant. 

During  the  development  of  the  constitutive  material  constants  it  was  foxmd  that  one 
set  of  variables  gave  an  excellent  match  of  the  initial  yield,  with  a  poor  prediction  of 
the  cyclic  response,  whilst  a  second  set  of  constants  grossly  misrepresented  the  initial 
yield,  but  gave  an  excellent  prediction  of  the  cyclic  response  (see  Figures  25  and  26). 
The  largest  difference  in  the  values  of  the  variables  between  the  two  scenarios  was 
shown  by  the  hardening  variables  fi  and  fi.  It  could  be  inferred  that  an  evolution  in  the 
value  of  the  variables  fi  and  within  the  model  between  the  initial  yield  and  the  cyclic 
condition  could  improve  the  overall  numerical  prediction.  In  effect  the  changes  to  the 
hardening  variables  would  not  be  evolutionary,  rather  a  one  off  shift  of  the  values  after 
the  material  has  initially  yielded.  This  would  allow  die  constitutive  model  to  give  a 
better  overall  prediction  of  the  material  behaviour,  as  it  would  accurately  model  the 
initial  yield  as  well  as  the  cyclic  condition. 

Another  material  characteristic  exhibited  by  D6ac  steel,  which  is  not  accurately 
predicted  by  the  constitutive  model,  is  the  relaxation  of  the  mean  stress  during  a  cyclic 
loading  application.  This  can  be  most  effectively  seen  for  a  CPLT  comparison  between 
the  experimental  results  and  the  numerical  prediction  as  in  Figure  24.  The  mean  stress 
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relaxation  is  a  direct  resrilt  of  the  maxitnum  compressive  stress  approaching  the 
maximum  tensile  stress  in  magnitude  as  the  material  progresses  cyclically.  During 
cyclic  loading  the  maximum  tensile  stress  remains  constant  whilst  the  maximum 
compressive  stress  relaxes,  thus  reducing  the  mean  stress  of  each  progressive  cycle. 
Most  unified  constitutive  models  can  accurately  predict  the  main  cyclic  development 
of  a  given  material,  however,  they  become  inadequate  when  it  comes  to  the  modeUing 
of  mean  stress  relaxation  under  strain  control,  or  conversely,  strain  shakedown  or 
ratcheting  under  stress  (load)  control  [14,15].  Shakedown  is  the  stabilisation  of  the 
mean  strain  imder  cyclic  loading,  whilst  ratcheting  is  the  divergent  behaviour  of  the 
mean  strain  as  the  number  of  cycles  increases. 

As  we  are  mainly  concerned  with  a  strain  controlled  loading  (CPLT),  mean  stress 
relaxation  is  of  interest  whilst  shakedown  and  ratcheting  are  not  analysed.  Basically, 
the  modelling  of  mean  stress  relaxation  is  a  direct  response  of  the  application  of 
kinematic  hardening  within  the  model.  A  linear  kinematic  hardening  rule  does  not 
predict  any  mean  stress  relaxation.  Whilst  the  application  of  a  non-linear  kinematic 
hardening  rule  does  predict  mean  stress  relaxation,  it  most  commonly  results  in  an 
over  prediction  [15].  As  can  be  seen  from  the  numerical  results  (Figure  24)  the  Stouffer 
constitutive  model  over  predicts  the  mean  stress  relaxation  which  is  a  common  trait  of 
such  a  constitutive  model.  Therefore,  it  is  unrealistic  to  expect  the  constitutive  model 
to  accxuately  predict  the  stress  relaxation  exhibited  by  D6ac  steel,  without  modification 
to  the  application  of  the  hardening  rule.  This  type  of  modification  would  require  re¬ 
formulation  of  the  constitutive  equations. 

It  should  be  noted  that  the  afore  mentioned  shortfalls  of  the  unified  constitutive  model 
apply  purely  to  the  material  behaviour  of  D6ac  steel.  The  tests  revealed  a  unique 
combination  of  material  properties  which  are  rmcommon  and  rmexpected  for  this  type 
of  high  strength  steel.  Consequently,  these  anomalies  cannot  be  effectively  predicted 
by  the  Stouffer  constitutive  model,  or  any  other  plasticity  modelling  technique.  Apart 
form  the  small  discrepancies  which  can  be  viewed  in  Figures  23  and  24,  the 
constitutive  model  predicts  the  cyclic  nature  of  the  D6ac  steel  well.  The  CPLT  loading 
cycle  is  of  a  non-symmetrical  nature  and  therefore  it  is  more  difficult  to  predict, 
especially  by  classical  plasticity  techniques.  So  although  the  unified  constitutive  model 
does  not  provide  an  exact  representation  of  the  behaviour  of  D6ac  steel,  it  does  predict 
the  cyclic  response  to  a  degree  which  is  acceptable  for  its  application  to  the  F-111 
DADTA  process. 


8.  Conclusions 


This  report  has  described  the  development  of  material  constants  of  D6ac  steel  for  the 
xmified  constitutive  model.  Test  data  was  obtained  through  a  serious  of  uniaxial  tests 
and  an  initial  prediction  of  the  material  constants  was  obtained.  These  initial 
predictions  were  refined  using  a  computer  representation  of  the  model  to  arrive  at  a 
final  set  of  values. 
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During  the  development  of  the  material  constants  a  number  of  material  anomalies 
were  observed  which  could  not  be  predicted  by  the  model.  These  anomalies  were 
peculiar  to  this  type  of  steel,  and  were  therefore  not  able  to  be  handled  by  the  unified 
model.  New  testing  techniques  were  also  developed  for  umaxial  specimens  which  led 
to  the  reduction  of  out  of  plane  bending  effects. 


Overall  a  set  of  material  constants  were  produced  which  allowed  the  cyclic  prediction 
of  the  material  behaviour  of  D6ac  steel.  These  results  are  to  be  applied  to  the 
calculation  of  the  residual  stress  field  in  the  F-111  WPF  undergoing  a  CPLT  load 


application. 
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Figure  3.  Tension/Compression  cycle  to  20,000  (jls 


Figure  4.  Compression/Fension  cycle  to  20,000  pLs 


Figure  5.  Tension/Compression  cycle  to  20,000  jus  at  10  fis/sec 
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Figure  6.  Tension/Compression  cycle  to  20,000  fisat  600  pis/sec 


Figure  7.  100%  CPLT  cycle 


Figure  8.  75 %  CPLT  cycle 
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Figure  12  ;  Specimen  3  Cycle  20  (60  microstrain/sec) 
Load  vs  Microstrain  (2  Clip  Gauges) 


Figure  12.  Specimen  3  Cycle  20  (60  microstrain/sec) 


Figure  13  ;  Specimen  4  (10  microstrain/sec) 
Load  vs  Microstrain 


Figure  13.  Specimen  4(10  microstrain/sec) 


Load  (kN) 
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Figure  15.  Specimen  6  (60  microstrain/sec) 


Load  (kN)  Load  (kN) 
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Figure  18  ;  Specimen  9  75%  CPLT  (60  microstrain/sec) 
Load  vs  Microstrain 


Figure  18.  Specimen  9  75%  CPLT  (60  microstrain/sec) 


Figure  19  ;  Specimen  1 0  50%  CPLT  (60  microstrain/sec) 
Load  vs  Strain 


Figure  19.  Specimen  10  50%  CPLT  (60  microstrain/sec) 
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STRESS  (MPa) 
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Figure  26.  Comparison  of  Numerical  and  Experimental  Response  for  Cyclic  Constants 
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Appendix  A: 

Uniaxial  Specimen  Finite  Element  Model 
Al.  Introduction 

The  specimens  used  to  calculate  the  unified  constitutive  state  variables  were  made  of 
D6AC  steel,  cylindrical,  with  a  length  (including  threads)  of  127mm  and  a  test  section 
diameter  of  9.5mm  (see  Figure  2).  All  of  the  tests  were  conducted  in  strain  control.  The 
controlling  extensometer  (25mm  gauge  length)  was  attached  across  the  test  section 
which  was  27mm  in  length.  A  fillet  curve  exists  at  the  junction  between  the  working 
section  and  the  threaded  region.  Given  that  the  overall  test  section  length  and  the 
gauge  length  are  quite  close,  and  that  the  knife  edges  of  the  extensometer  sit  extremely 
close  to  the  fillet,  a  stress  analysis  of  this  region  was  performed  to  investigate  any 
possible  interference  to  the  extensometer. 

To  analyse  the  stress  distribution  within  the  specimen  a  detailed  finite  element  model 
was  constructed.  The  model  consisted  of  the  specimen  only  and  was  aimed  at 
providing  details  of  the  transfer  of  load  from  the  threaded  connections  through  to  the 
test  section.  In  addition,  the  through  thickness  stress  distribution  of  the  specimen  was 
examined.  It  was  anticipated  that  high  regions  of  stress  would  appear  at  the  section 
change  between  the  threaded  area  and  the  test  section,  whilst  the  stress  in  the  test 
section  would  be  constant. 

The  initial  distortion  of  the  specimen  due  to  machining  and  heat  treatment  has  been 
ignored  in  the  analysis. 


A2.  Finite  Element  Model  Description 

The  finite  element  model  of  the  uniaxial  specimen  was  produced  using  the  PAFEC 
Finite  Element  Package  currently  in  use  at  AMRL.  The  model  contained  4483  nodes, 
720  20-noded  isoparametric  brick  elements,  360  15-noded  isoparametric  prism 
elements  and  13,041  degrees  of  freedom.  A  uniform  mesh  was  employed  with 
increased  detailed  meshing  throughout  the  working  section  of  the  specimen  (see 
Figure  Al).  The  model  was  restrained  at  one  end  and  loaded  at  the  other,  thus 
simulating  the  way  in  which  the  specimen  was  loaded  in  the  testing  machine.  To 
simulate  the  correct  level  of  restraint,  die  nodes  around  the  periphery  of  the  threaded 
area  were  restrained  in  all  directions  thus  emulating  a  clamped  end.  The  loading  of  the 
opposite  end  was  applied  as  a  distributed  tensile  load. 

The  general  application  of  load  is  performed  by  varying  the  distribution  along  the 
threaded  end  whereby  most  of  the  load  is  applied  to  the  first  few  threads  and  then 
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rapidly  diminishes  along  the  remaining  length  of  the  thread.  For  this  model  a 
mathematical  estimate  of  this  distribution  was  employed  which  was  based  around  the 
following  equation: 


P(z) 


Pm 

sinh(m//) 


cosh(/Mz) 


(Al) 


where  p(z)  is  the  force  taken  up  by  the  z*'’  turn  of  the  thread,  P  is  the  total  load  acting 
along  the  thread,  m  is  a  coefficient  dependant  upon  the  construction  of  the  thread,  and 
H  is  the  length  of  engagement  [16]. 

The  model  was  loaded  with  a  equivalent  distributed  axial  tensile  load  of  llOkN. 
Although  the  specimen  experiences  plasticity  at  this  level  of  loading,  a  plastic  analysis 
was  not  performed.  The  loading  was  applied  as  point  loads  arotmd  the  periphery  of 
the  threaded  section  of  the  specimen,  with  the  longitudinal  distribution  conforming 
with  that  previously  discussed. 


A3.  Results 

As  expected  the  highest  regions  of  stress  were  found  to  exist  in  the  fillet  radius  where 
the  two  sections  meet,  occurring  at  the  beginning  of  the  working  section.  The  stress 
plots  (Figure  Al)  depict  the  surface  stress  distribution  throughout  the  specimen,  and 
highlight  the  maximum  stress  occurring  at  the  radial  section  change.  Figure  A2  is  a 
plot  of  stress  in  the  longitudinal  direction.  The  resulting  elastic  stress  concentration 
factor  at  the  surface  of  the  specimen  was  calculated  to  be  1.2.  Figure  A3  is  a  plot  of  a 
section  cut  through  the  middle  of  the  specimen.  It  can  be  seen  in  this  plot  that  the  areas 
of  maximum  stress  exist  around  the  periphery  of  the  specimen  and  do  not  propagate 
through  the  thickness.  From  this  plot  it  is  also  interesting  to  note  that  the  specimen 
achieves  a  constant  through  the  thickness  loading  throughout  the  length  of  the  test 
section. 


A4.  Conclusion 

The  results  of  this  analysis  clearly  show  that  the  uniaxial  specimen  experiences  a 
constant  through  the  thickness  loading  throughout  the  working  section.  In  addition,  it 
can  be  seen  that  localised  areas  of  high  stress  exist  at  the  section  change.  As  this  is 
where  the  clip  gauge  knife  edges  sit,  it  is  possible  that  the  clip  gauge  may  be  reading 
slightly  inconsistent  results.  However,  since  the  localised  stress  intensity  factor  is  1.2, 
and  the  stress  distribution  is  uniform  across  the  working  section,  it  is  anticipated  that 
the  effect  on  the  experimental  data  would  be  negligible. 
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Figure  Al.  Uniaxial  Finite  Element  Mesh 


Figure  A2.  Tension  Load  (llOkN)  Stress  Distribution 


LINEAR 


Figure  A3.  Tension  Load  (llOkN)  Stress  Distribution  X-Section 
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